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Measurement of the forward-backward asymmetry in top-antitop quark 
production in proton-antiproton collisions 

D. Orbaker for the D0 collaboration 
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We present a new measurement of the forward-backward asymmetry in tt production in pp col- 
lisions in the i-|-jets channel at DO. We measure asymmetries from two different observables and 
unfold the data to allow the results to be compared to standard model predictions. For the unfolded 
T— I asymmetry based on the Aj/ variable, we measure (19.6 ± 6.5)%, compared with an MC@NLO-based 

^""l prediction of (5.0±0.1)%. We also discuss the correlation between the asymmetry and the transverse 

momentum of the tt system. 

(N 



INTRODUCTION 



o 
O 

In this note we present a new measurement of the forward-backward asymmetry (Ape) in ti production in 
^""^ pp collisions using the /-|-jets channel at DO [T]. Using 5.4 fb^^ of data, this measurement is an update to the 
first DO result using 0.9 fb~^ |2]. Recently, the CDF collaboration released a result presenting an asymmetry 
^ that is a 2-3 sigma deviation from the inclusive standard model (SM) prediction and has a dependence on the 
O invariant mass of the tt system, rutt [3] • 

^ In the SM, the main contribution to A-pB in tt events comes from quantum chromodynamics (QCD). The first 
calculation for asymmetry at the Tevatron, made more than a decade ago, predicts an A-pB of ^ 5% [4]. More 
<~| recent calculations, some of which include electroweak (EW) effects, predict asymmetries up to ^ 9% [5H7]. 
Recent theoretical efforts have been made to include the asymmetry in a framework of measurements including 
I top quark polarization and spin correlation [8] . 
^ To compare our results with theory we present two different types of results. The first type are reconstruction 

fS| level results. Reconstruction level asymmetries appear in the DO detector, after the effects of detector acceptance 
and event reconstruction. Quantities measured at the reconstruction level cannot be directly compared with 
theory. For this reason we also present the measurement at the production level, before effects of selection and 
reconstruction. To access the production level for a given quantity we unfold the distributions in data after 
background subtraction. 

^ II. THE DO DETECTOR 

> 

k> DO is a general purpose particle detector [S] ■ Going outward from the beam pipe, the first subsystem of the 

detector is the central-tracking system, which consists of a silicon microstrip tracker (SMT), as well as a central 
C3 fiber tracker (CFT). Both portions of the tracker are located inside a 2 T superconducting solenoid magnet. 
The designs are optimized to track and vertex charged particles at rapidities up to jr;! < 3 for the SMT and 
jryl < 2.5 for the CFT. Right outside of the solenoid are central and forward preshower detectors. Three liquid 
argon and uranium calorimeters sit outside of the preshower: one central section (CC), which covers |?7| up to 
ss 1.1 and two end sections (EC) that extend the range in |?7| to ~ 4.2. All three sections of the calorimeter are 
housed in separate cryostats [TD]. Furthest away from the collision region is the muon system, which covers a 
range of I77I < 2 and consists of a layer of tracking detectors and scintillator trigger counters, followed by 1.8 T 
toroids and then two similar detector and scintillator layers [H]. Plastic scintillator arrays placed in front of 
the EC cryostats measure the luminosity. Both the trigger and data acquisition systems are designed to handle 
the high instantaneous luminosities of the Tevatron Run II. 



III. EVENT SELECTION AND RECONSTRUCTION 



The requirements for events to be selected are very similar to those found in the most recent DO ^+jets 
crosssection measurement [l^]. In the /-|-jets channel, the top quark pair decays as such: ti W^bW^b, 
with one W boson decaying leptonically, W ^ Iv and the other W boson decaying hadronically W ^ qq . 
Looking at ti bblvqq , one can see that there are six decay products in the final state: a lepton, which in this 
measurement is either a muon (/i) or an electron (e), missing transverse energy from the neutrino {fr^) and four 
jets, two of them jets originating from b quarks. 
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Keeping the six decay products in mind, we require four jets with transverse momentum {j)t) greater than 
20 GeV and pseudorapidity {\rj\) less than 2.5. At least one of these jets must pt > 40GeV. We require one 
and only one lepton; either an electron with pt > 20GeV and |r/| < 1.1 or a muon with pj- > 20GeV and 
\ri\ < 2.0. Missing transverse energy must be greater than 20 GeV in the e+jets channel and 25 GeV in the 
/Lt+jets channel. 

Jets originating from b quarks are identified by displaced vertices via a 6-tagging neural net algorithm [13j . 
Every event must have at least one jet that passes the ^-tagging requirements. The efficiency for a b jet to pass 
the 6-tagging requirements is ^70%, while likelihood of a jet originating from quarks of lighter weight to mimic 
a 6-jet and be misidentified is ^8%. For interested readers, more detail may be found here [H]. 

We employ a constrained kinematic fit to reconstruct the six decay objects into the top quark four vectors 
under the tt decay hypothesis [TS]. A function based on the detector resolution is minimized for each jet- 
parton permutation and only the assignment with the lowest is kept. To reduce the number of possible 
assignments, at least one of the jets associated with a b quark by the kinematic fit must be 6-tagged. A total of 
four constraints are used: each reconstructed top quark must have a mass of 172.5 GeV and each reconstructed 
W boson must have a mass of 80.4 GeV. The jet-parton assignment is correct ~ 70% of the time. 



IV. SIMULATION AND BACKGROUNDS 



Events from ti decays are simulated with mc@nlo, with herwig used for showering (TBI [I?]. Background 
events from W boson produced in association with jets (VF+jets) decays are simulated in a similar fashion 
with ALPGEN+PYTHIA [THl US] • We also take into account background events from multijet (MJ) production, 
where one jet mimics the signature of a lepton. Other backgrounds such as single top, diboson and Z+jets 
are insignificant and not considered in this analysis. To be comparable to data, events from the Monte Carlo 
simulations are run through the DO detector simulation '20] and the reconstruction sequence used on data. To 
evaluate the contribution from the MJ background we use control samples from DO data. 



V. DEFINITION OF ASYMMETRIES 

The forward-backward asymmetry is one of the quantities useful for summarizing a differential distribution. 
The forward-backward asymmetry is the difference between the fraction of events defined as forward {Nf) and 
the fraction of events defined as backward {Nb), 

APB = (1) 

We present the asymmetry in two different observables; Ay = yt — Vt = <li{yt,iep ~ yt,had) and qiyi, where 
the rapidity is defined as y = ^In (^ ^-p' ) ' l^- charge of the lepton, t, lep is the leptonically decaying top 

quark and t, had is the hadronically decaying top quark. 

For Ay, the asymmetry is 

^ A^(Ay > 0) - jV(Ay < 0) 
™ iV(Ay > 0) + 7V(Ay < 0) ■ ^' 

The forward-backward asymmetry for qiyi, the so-called lepton-based asymmetry, is 

j^l ^ Njqm > 0) - Njqiyi < 0) 
™ Niqiyi>0) + N{qiyi<0)- ^' 



VI. PREDICTED ASYMMETRIES 



We predict Aps and Apg at the reconstruction level using the events generated from MC@NLO and run 
through the DO detector simulation with the same selection criteria and reconstruction used on data. The 
predicted asymmetries for both the reconstruction and production levels are summarized in Table [T] Production 
level asymmetries are taken before events are run through through the DO detector simulation. 
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TABLE I: Predictions from the MCQNLO event generator. 





Channel Ays (%) 


A^B (%) 


Generated 
Reconstructed 


Z+jets 5.0 ±0.1 

e+jets 2.4 ±0.7 
/x±jets 2.5 ± 0.9 
/±jets 2.4 ±0.7 


2.1 ±0.1 

0.7 ±0.6 
1.0 ±0.8 
0.8 ±0.6 



VII. RECONSTRUCTION LEVEL TECHNIQUE AND RESULTS 



Both the amount of events from signal and background processes, i.e. the sample composition, and the 
asymmetry from ti events are measured via a template-based maximum likelihood (ML) fit. A discriminant 
made up of four input variables shown in Figure [T] that are loosely correlated with Ay is trained to separate the 
W+jets background from the ti signal. Four different templates containing the discriminant and the sign of the 
asymmetry are fit to the data: a signal template made up with events where Ay > 0; another signal template 
made up with events where Ay < 0; the template for the W+jets background with the asymmetry taken from 
simulation; the template for the MJ background with the discriminant shape and asymmetry taken from the 
control data sample. 

The results using the method described in the previous paragraph for asymmetries from both Ay and qiyi are 
presented in Tables |TT] and [lllj respectively. For the lepton-based measurement, an additional selection criterion 
of \yi\ < 1.5 is applied. To search for new physics, the asymmetry is also presented for different invariant mass 
regions and different magnitudes of Ay in Table llV 





FIG. 1: Inputs to the discriminant used for the reconstruction level measurement, a) from the kinematic fit. b) 
Transverse momentum of the fe-tagged jet with the largest transverse momentum, c) The minimum relative momentum 
between any two of the hardest four jets, d) Invariant mass of the two jets assigned by the kinematic fit as daughters of 
the hadronically decaying W boson. 



VIII. UNFOLDING TECHNIQUE AND RESULTS 



To infer the asymmetry at the production level, we use regularized unfolding to correct the reconstruction 
level Ay distribution, which is the Ay distribution in data after subtracting the contribution from the back- 
grounds. Unfolding corrects for effects from detector reconstruction and acceptance, which change Ay at the 
reconstruction level. In order to take into account the fine details of the reconstruction the unfolding procedure. 
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Discriminant Discriminant 

FIG. 2: The distribution in the discriminant for data and simulation for events with Ay < (a) and Ay > (b). 

TABLE II: Numbers of events in data, results of fits for sample composition and ^fb, and predictions for j4fb. The 
asymmetries are given at reconstruction level, with their total uncertainties. The sample compositions are given with 
the fit uncertainties. 





l+>4: jets 


e+>4 jets 


/i+>4 jets 


1+4 jets 


l+>5 jets 


Raw Nf 


849 


455 


394 


717 


132 


Raw Nb 


732 


397 


335 


597 


135 




1126±39 


622±28 


502±28 


902±36 


218±16 


Nw+}cts 


376±39 


173±28 


219±27 


346±36 


35±16 




79±5 


56±3 


8±2 


66±4 


13±2 




9.2±3.7 


8.9±5.0 


9.1±5.8 


12.2±4.3 


-3.0±7.9 


MC@NLO Afb (%) 


2.4±0.7 


2.4±0.7 


2.5±0.9 


3.9±0.8 


-2.9±1.1 



based on a modified version of ROOT's TUnfold class, uses a migration matrix with 50 bins for the reconstruc- 
tion level and 26 bins for the production level 21 . The acceptance correction multiples each bin by -^rjj-: where 
effi is the selection efficiency for bin i in Ay. Regularization smoothes the large statistical ffuctuations between 
bins in the Ay distribution. The unfolded Ay distribution is summarized in the two bins of the asymmetry. 

In addition to regularized unfolding, we also perform a four-bin maximum likelihood unfolding of the Ay 
distribution. For the four-bin procedure, the boundaries of the bins are at Ay = —3, —1, 0, 1, 3. The result of 
the four-bin unfolding procedure is Afb = (16.9 ± 8.1)%. 

Because migrations in between bins for the lepton-based are very small, only an acceptance correction is used 
to unfold the qiyi distribution. Results from the unfolding are shown in Table |Vj 



TABLE III: Numbers of events in data, results of fits for sample composition and Ap-g, and predictions for Ap^- The 
asymmetries are given at reconstruction level, with their total uncertainties. The sample compositions are given with 
the fit uncertainties. 





l+>4: jets 


e±>4 jets 


/x±>4 jets 


/±4 jets 


/±>5 jets 


Raw Np 




867 


485 


382 


730 


137 


Raw Ng 




665 


367 


298 


546 


119 


Nu 




1096 ± 39 


622 ± 28 


474 ± 27 


881 ± 36 


211 ± 16 


Nw+jets 




356 ± 39 


173 ± 28 


198 ± 27 


323 ± 36 


31 ± 16 


Nm.1 




79 ±5 


56 ±3 


8±2 


66 ±4 


14 ± 2 


A'pB (%) 




14.2 ±3.8 


16.5 ±4.9 


9.8 ±5.9 


15.9 ±4.3 


7.0 ±8.0 


MC@NLO ApB 


(%) 


0.8 ±0.6 


0.7 ±0.6 


1.0 ±0.8 


2.1 ± 0.6 


-3.8 ± 1.2 
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TABLE IV: Reconstruction-level Afb by subsample. 





ApB (%) 


Subsample 


Data 


MC@NLO 


niu < 450 GeV 
mtt > 450 GeV 
|Ay| < 1.0 
lAyl > 1.0 


7.8 ±4.8 
11.5 ±6.0 
6.1 ±4.1 
21.3 ±9.7 


1.3 ±0.6 

4.3 ± 1.3 

1.4 ±0.6 
6.3 ± 1.6 




FIG. 3: Distribution in qiyi (a) and Ay (b) for data and simulation. The sample composition is taken from the fit. 

IX. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties, shown in Tables |VI| and |VIH are computed for both the reconstruction and 
production level measurements, as well as for the reconstruction level predictions. Note that although we make 
a detailed accounting of the systematics effects, the measurement is dominated by statistical uncertainty. The 
uncertainties into seven different categories. Further details about each category of uncertainty may be found 
in Ref [1]. 

X. CROSS CHECKS 

We perform multiple cross checks to ensure that the measurement is accurate. The asymmetry from the 
W-|-jets background is taken from simulation, specifically ALPGEN±PYTHIA. To check that the simulated asym- 
metry is accurate, we use a similar template-based measurement which includes events without a 6-tagged jet, 
as seen in Figure [Sj but meeting every other event selection criteria. These events are dominated by the W±jets 
background. The fitted asymmetry for W±jets using this method is in good agreement with the simulated 
asymmetry. 

The polarities of the solenoid and toroid magnets of the DO detector are regularly flipped. To ensure that there 
is no inherent bias in the asymmetry from the detector, measurements are made for each of the four polarity 
settings for both the Ay- and lepton-based asymmetries. No significant differences were found between samples 
with different settings of magnet polarities. Similarly, breaking the lepton-based asymmetry into samples with 
positive and negative charge did not significantly affect the result. 



TABLE V: Ay- and lepton-based Asymmetries. 





Afb (%) 


^FB (%) 




Reconstruction level 


Production level 


Reconstruction level 


Production level 


Data 


9.2±3.7 


19.6±6.5 


14.2 ± 3.8 


15.2 ± 4.0 


MC@NLO 


2.4±0.7 


5.0±0.1 


0.8 ± 0.6 


2.1 ± 0.1 
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Produced Ay 

FIG. 4: The 50 by 26 bin migration matrix used for the regularized unfolding. 



TABLE VI: Absolute uncertainty on ylpB (%) 





Reco. level 


Prod, level 


Source 


Prediction 


Measurement 


Measurement 


Jet reco 


±0.3 


±0.5 


+ 1.0 


JES/JER 


+0.5 


-0.5 


-1.3 


Signal modeling 


±0.3 


±0.5 


+0.3/-1.6 


6-tagging 




±0.1 


+0.1 


Charge ID 




+0.1 


+0.2/-0.1 


Bg subtraction 




±0.1 


+0.8/-0.7 


Unfolding Bias 






+1.1/-1.0 


Total 


±0.7/-0.5 


+0.8/-0.9 


+1.8/-2.6 



XI. DISCUSSION 

In summary, the reconstruction level results for asymmetries based on Ay and qiyi are presented in Tables [Tl| 
and |III[ respectively. Results from unfolding to compare to the production level predictions are shown in 
Table |V| The measured values for both asymmetries are higher than the predictions MC@NLO. 

In addition to the measurements, we investigated the behavior of asymmetry on gluon radiation via the 
transverse momentum of the tt system, p". In the SM, the asymmetry arises at next-to- leading-order from an 
interference of various diagrams. Diagrams which contain external gluon lines, either in the initial state or the 
final state, interfere to produce a negative asymmetry. On the other hand, the interference between diagrams 
in which all gluon lines are internal produces a positive asymmetry. 

Because is made up of the six detected objects from tt decay, the recoil from low energy gluon radiation 
can be measured. The dependence of the asymmetry across the spectra contains both soft and hard gluon 
radiation, and the behavior is shown in fine detail in Figure tel The dependence of Afb on simulated 



TABLE VII: Absolute uncertainty on Ape (%) 





Reco. level 


Prod, level 


Source 


Prediction 


Measurement 


Measurement 


Jet reco 


+0.3 


+0.1 


+0.8 


JES/JER 


+0.1 


-0.4 


+0.1/-0.6 


Signal modeling 


+0.3 


+0.5 


+0.2/-0.6 


f)-tagging 




+0.1 


+0.1 


Charge ID 




+0.1 


+0.2/-0.0 


Bg subtraction 




+0.3 


+0.6 


Total 


+0.5 


+0.7 


+1.0/-1.3 
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■ ft 

□ IV+jets 



■ Multijet i 
• 1 4/ \ t 



D0, 5.4 fb"^ 
(a) 





FIG. 5: Distributions in qiyi (a) and Ay (b) for events with zero ^-tagged jets. The sample composition is taken from 
the fit. 



with MC@NLO and pythia is shown. Two different dependences are shown for pythia, one with the angular 
coherence (angular ordering) parameter on (shown for the curve corresponding to PYTHIA 6.425 D6-Pro). The 
other curve, for the pythia 6.425 SOA-Pro setting, has angular coherence turned off, and the dependence of 
A-pB on is no longer present. We note that we only use pythia for qualitative purposes, and mc@nlo is 
used for all predictions. Because the selection efficiency versus is not constant (Figure [H]), we turn off the 
dependence of A-p-Q and for one of the systematics. 

We briefly mention the comparison of data and simulation for p^. In Figure 7| we show the distribution in 
for both MC@NLO and pythia with initial state radiation (ISR) turned off. Here pythia with modified initial 
state radiation does a better job describing the data. Two things should be noted. One, the resolution for p" 
is poor. Bins widths are shown with half the resolution, meaning two bins is about the best it can get. Two, 
the distribution in the number of jets, which agrees very well with data using mc@nlo, but no longer agrees 
while when using pythia with ISR turned off. We present no concrete conclusion for this issue, but show what 
we have found so far. 



"j: 0.3 
0.2 
0.1 

OF- 
-0.1^ 
-0.2^ 
■0.3^ 
-0.4^ 




MC@NLO 3.4 
PYTHIA 6.425 SOA-Pro 
PYTHIA 6.425 D6-Pro 



10 20 



30 _ 40 50 60 70 80 90 
tt transverse momentum [GeV] 




10 20 30 40 50 60 70 80 90 100 

tt transverse momentum [GeV] 



FIG. 6: ApB versus for mcQnlo (black) and pythia with different settings (a). The selection efficiency versus p" 
(b). 



D0, 5.4 fb'^ tt MC@NLO 

(a) □ IV+jets 

Multijet 
• Data 




D0, 5.4 fb"' tt PYTHIA ISR off 

(b) □ W+jets 

Multijet 
• Data 



10 20 30 40 _ 50 60 70 80 90 100 

Reconstructed ft transverse momentum [GeV] 




10 20 30 40 _ 50 60 70 80 90 100 

Reconstructed tt transverse momentum [GeV] 



FIG. 7: Distribution in p" for data and simulation for mcOnlo (a) and pythia with initial state radiation turned off 
(b). 



8 



Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 



XII. CONCLUSION 



We present measurements of the forward-backward asymmetry in ti events using both Ay and qiyi. At the 
reconstruction level, the Ay-based asymmetry is (9.2 ± 3.7)% in data and (2.4 ± 0.7)% for the MC@NLO-based 
prediction. The lepton-based reconstruction level asymmetry is (14.2 ± 3.8)% for data and (0.8 ± 0.6)% for 
the MC@NLO-based prediction. 

In addition to the reconstruction level measurements, we also present unfolded asymmetries for Ay and qiyi- 
Unfolding the Ay distribution gives an asymmetry of (19.6 ± 6.5)%, compared to MC@nlo production level 
predictions of (5.0 ± 0.1)%. In comparison, unfolding the qiyi distribution results in an asymmetry of (15.2 ± 
4.0)%, to be compared with (2.1 ± 0.1)%. We note that other SM predictions measure larger asymmetries. 

The measurements from DO data are consistently higher than the predictions. We mentioned some potential 
limitations of these predictions. 
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